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ABSTRACT 

In Paper I (Fiege & Pudritz, 1999), we constructed models of filamentary molecular 
clouds that are truncated by a realistic external pressure and contain a rather gen- 
eral helical magnetic field. We address the stability of our models to gravitational 
fragmentation and axisymmetric MHD-driven instabilities. By calculating the domi- 
nant modes of axisymmetric instability, we determine the dominant length scales and 
growth rates for fragmentation. We find that the role of pressure truncation is to de- 
crease the growth rate of gravitational instabilities by decreasing the self-gravitating 
mass per unit length. Purely poloidal and toroidal fields also help to stabilize fila- 
mentary clouds against fragmentation. The overall effect of helical fields is to stabilize 
gravity- driven modes, so that the growth rates are significantly reduced below what 
is expected for unmagnetized clouds. However, MHD "sausage" instabilities are trig- 
gered in models whose toroidal flux to mass ratio exceeds the poloidal flux to mass 
ratio by more than a factor of ^ 2. We find that observed filaments appear to lie 
in a physical regime where the growth rates of both gravitational fragmentation and 
axisymmetric MHD-driven modes are at a minimum. 

Key words: molecular clouds - MHD - instabilities - magnetic fields. 
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1 INTRODUCTION 

It has long been known that filamentary molecular clouds 
often have nearly periodic density enhancements along 
their lengths. One has only to look at the Schneider and 
Elmegreen (1979) catalogue of filaments to observe many ex- 
amples of so-called globular filaments, where quasi-periodic 
cores are apparent. Such periodicity has also been noted 
along the central ridge of the T-shaped filament of Orion A 
(Dutrey et al. 1991). The fragmentation of self-gravitating 
filaments has accordingly received some attention. Given 
that overdense fragments may eventually evolve into star- 
forming cores, it is important to have as complete a picture 
of fragmentation as possible. 

Chandrasekhar and Fermi (1953) first demonstrated 
that uniform, incompressible filaments threaded by a purely 
poloidal magnetic field are subject to gravitational instabil- 
ities, and that magnetic field helps to decrease the growth 
rate of the instability. More recently, other authors (Naga- 
sawa 1987, Nakamura, Hanawa, & Nakano 1993, Tomisaka 
1996, Gehman 1996, etc.) have studied the stability of com- 
pressible, magnetized filaments. We perform a stability anal- 
ysis for the equilibrium models of filamentary molecular 



clouds that we described in Fiege & Pudritz 1999 (here- 
after FPl). Our equilibria are threaded by helical magnetic 
fields and are truncated by the pressure of the ISM. Most im- 
portantly, our models are constrained by observational data 
and are characterized by approximately r~^ density profiles, 
which are in good agreement with the observational data 
(Alves et al. 1998, Lada, Alves, and Lada 1998). 

We demonstrated in FPl that most filaments are well 
below the critical mass per unit length required for radial 
collapse to a spindle, and that any magnetized filamentary 
cloud that is initially in a state of radial equilibrium is stable 
against purely radial perturbations. Here, we consider more 
general modes of axisymmetric instability to find the growth 
rates and length scales on which our models may break up 
into periodic fragments along the axis. We show that our 
models are subject to two distinct types of instability. The 
first is a gravity-driven mode of fragmentation, which occurs 
when the toroidal field component is relatively weak com- 
pared to the poloidal. The second type of instability arises 
when the toroidal magnetic field is relatively strong com- 
pared to the poloidal field, and probably represents the ax- 
isymmetric "sausage" instability of plasma physics. We find 
a stability criterion for these modes, which we discuss in Sec- 
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tion 0. We find that the gravitational and MHD modes blend 
together at wavelengths that are intermediate between the 
long wavelengths of purely gravitational modes, and the very 
short wavelengths of purely MHD-driven instabilities. These 
intermediate modes, in fact, represent a physical regime in 
which filamentary clouds fragment very slowly. 

Once the process of fragmentation begins, it becomes 
non-linear on a timescale governed by the growth rate of 
the instability. We stress that our analysis is linear and, 
therefore, can say nothing about the non-linear evolution 
of filamentary clouds. However, the linear stability is impor- 
tant because it determines the dominant length scale for the 
separation of fragments, and the largest scale for substruc- 
ture in filamentary clouds. These fragments might eventually 
evolve into cores In the present work, we limit our discus- 
sion to how the magnetic field configuration and external 
pressure affect spacing and mass scale of fragments. 

We consider only axisymmetric modes here. The analy- 
sis of non-axisymmetric modes, most notably the kink insta- 
bility (cf. Jackson 1975), will be investigated in a separate 
paper. One might expect fast-growing kink instabilities be- 
cause of the toroidal character of the outer magnetic field. 
However, the centrally concentrated poloidal field is much 
stronger than the peak toroidal field in most of our mod- 
els; this magnetic "backbone" should largely stabilize our 
models against the kink instability. 

A brief plan of our paper is as follows. In Section bl 
we discuss the equilibrium state of the molecular filament 
and atomic envelope, as well the general method of our sta- 
bility calculations. We derive the equations of motion and 
all boundary conditions in Sections and 0. We discuss the 
general stability properties of filaments in Sections, and the 
test problems that we used to verify our numerical code in 
Section g. We discuss the fragmentation of purely hydro- 
dynamic filaments in Section W, and filaments with purely 
poloidal or toroidal magnetic fields in Section H. In Section bl 
we analyze the stability of the observationally constrained 
models with helical fields discussed in Paper I. Finally, in 
Section lid, we discuss the significance of our results, and 
summarize our main findings. 



2 GENERAL FORMULATION OF THE 
PROBLEM 

The models described in FPl involve three parameters, two 
to describe the mass loading of the magnetic flux lines, and 
a third to specify the radial concentration of the filament. 
A Monte Carlo exploration of our parameter space led us 
to construct a class of magnetostatic models that are con- 
sistent with the observations. We use the equations of lin- 
earized MHD to superimpose infinitesimal perturbations on 
these observationally allowed equilibria. In general, we as- 
sume that molecular filaments are embedded in a less dense 
envelope of HI gas. In addition, we assume that the envelope 
is non-self-gravitating and is in equilibrium with the grav- 
itational field of the filament. Thus, the envelope is most 
dense at the interface between the molecular and atomic 
gas, and slowly becomes more rarified with radial distance. 
There is, in fact, at least one known example of a filament 
that is embedded in molecular gas of lower density. The f- 
shaped filament of Orion A is a very dense filament that is 



presently fragmenting and forming stars. Since its external 
medium is molecular, our treatment of the external medium 
as non-self-gravitating strictly does not apply to this case. 
We also assume that the poloidal field component (Bz) 
remains constant in the HI envelope, while the toroidal field 
component (-B^) decays with radius as r"^, so that the equi- 
librium state of the atomic gas is current free. We note that 
the magnetic field exerts no force in zeroth order in this con- 
figuration; thus, equilibrium is determined by hydrostatic 
balance alone. We model t he H I envelope as a polytrope. 



which we discuss in Section 2.1 



Of the authors mentioned in section hi, only Chan- 
drasekhar and Fermi (1953) and Nagasawa (1987) consid- 
ered filaments that are truncated at finite radius. Chan- 
drasekhar and Fermi considered the external medium to 
be a vacuum, while Nagasawa considered an infinitely hot 
and non-conducting external medium of zero density but fi- 
nite pressure. We treat the external medium as a perfectly 
conducting medium of finite density. By using the most 
general possible boundary conditions between two perfectly 
conducting media, we self-consistently solve the equations 
of linearized MHD in both the molecular filament and the 
surrounding HI envelope. Our approach allows us to both 
study the effects of a finite density envelope on the predicted 
growth rates of the instability, and also to predict what gas 
motions might arise in the envelope during fragmentation. 



2.1 The Polytropic HI Envelope 

The analysis of Paper I provides the equilibrium state of the 
self-gravitating molecular filament. However, we must also 
find an equilibrium solution for the HI envelope. We assume 
that the envelope is non-self-gravitating and current-free, 
with both field components continuous across the interface. 
Thus, the field in the envelope is given by 



B. 
B</) 



= B. 



B, 



<t>s 



(^\ 



(1) 



\RsJ ' 

where Rs is the radius of the molecular filament, as defined 
in Paper I, and B^s and B^s are the field components at 
the surface. 

It is required that the the total pressure of the gas and 
magnetic field balance at the interface. Thus the boundary 
condition on the total pressure is given by 



P-f- 



Bi + Bi 



Stt 



0, 



(2) 



where we use square brackets to denote the jump in a quan- 
tity across the interface. Since we have assumed that the 
field components are continuous in our model, the pressure 
must be continuous also. 

The magnetic field exerts no force on the gas in our 
current-free configuration. Therefore, we may derive the 
equilibrium structure of the envelope by taking the gas to 
be in purely hydrostatic balance with the gravitational field 
of the molecular filament; 



= 



dr 






(3) 



where subscript e refers to the HI envelope, and all quanti- 
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ties have their usual meanings. We assume that the gas is 
polytropic with polytropic index 7e; 



Ps 



_Pe_ 
Pe,S 



(4) 



The gravitational potential outside of a cylindrical filament, 
with mass per unit length m, is given by 



^ = 2Gm\n 



Rs' 



(5) 



where we choose the zero point of $ to be at the surface of 
the filament (radius Rs). With the help of kl and H, it is easy 
to solve for the density structure of the HI envelope: 



Ps 



o-Ite 



L r \'"' 

-GTTiln -— — I 



(6) 



The pressure distribution may then be obtained from equa- 
tion W. We find that the pressure falls very slowly with ra- 
dius, especially for 7e < 1. Generally we will take 0.25 < 
7e < 0.5 to simulate an atomic envelope that is nearly iso- 
baric out to very large radius. 



3 THE EQUATIONS OF LINEARIZED MHD 

In this Section, we set up the eigensystem of equations gov- 
erning the motion of a perfectly conducting, self-gravitating, 
magnetized gas. We assume that all perturbed variables can 
be expressed as Fourier modes. For example, we write the 
perturbed density pi in the form 

pi(r,z,t) = pi(r)e'("*+™^+''-^\ (7) 

where m — for the axisymmetric modes considered in this 
paper. We note that all quantities are written in dimension- 
less form using the dimensional scalings introduced in Paper 
I. 

We shall customarily reserve the subscript "0" for equi- 
librium quantities and "1" quantities in first order perturba- 
tion. In Appendix IA^, we show that the linearized equations 
of MHD reduce to the following set of coupled differential 
equations. The momentum equation combined with induc- 
tion and continuity equations becomes 



UJ povi = 7V • (povi) 



+ 



+ 



-<i(VxBo)x 



Vx 



VxVx 



Bo 

PoViX 

Po 



Bo 

POVlX 

Po 



xBo 



- poV- (itj$i) -I- V$oV- (povi), 
and Poisson's equation becomes 

= vVi + e^v- (povi). 



(8) 



(9) 



where tpi = iuj^i. The adiabatic index 7 is allowed to be 
discontinuous at the interface between the molecular and 
atomic gas; in general we define 



J: 



1, r<Rs 
7e, r > Rs 



(10) 



We have also included a parameter Q in equation M that 
allows us to "turn off" self-gravity in the HI envelope, as 
discussed in section H: 
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g = 



1, r< Rs 
0, r>Rs 



(11) 



Equations M and H represent the eigensystem, written 
entirely in terms of the perturbed momentum density povi 
and <pi . Appendix M shows that equations H and M can be fi- 
nite differenced and written in the form of a standard matrix 
eigenvalue problem: 



L*, 



(12) 



where L is a 4 x 4 block matrix operator, and ^ is an "eigen- 
vector" made up of the components of the momentum den- 
sity and the modified gravitational potential <p: 



* 



PO«rl 
POfcj,! 
pOVzl 

<Pl 



(13) 



It has been shown by Nakamura (1991) that the equations 
of self-gravitating, compressible, ideal MHD are self-adjoint; 
thus the eigenvalue —cu^ of equation [12 is purely real, and 
Lo may be purely real or purely imaginary. When lu is real, 
the perturbation is stable, since it just oscillates about the 
equilibrium state. On the other hand, when a; is imaginary, 
the perturbation grows exponentially, and the system is un- 
stable. It is these unstable modes that we are primarily con- 
cerned with, since they lead to fragmentation. 



4 BOUNDARY CONDITIONS 

Equation y_2| does not yet fully specify the eigensystem since 
boundary conditions have yet to be imposed. There are three 
radial boundaries in our problem; the first two are obviously 
the boundaries that occur at r = and r ~* 00. The third 
is the internal boundary that that separates the molecular 
filament and the HI envelope. As we discussed in Section bl 
we cannot directly difference across the interface because the 
density is discontinuous on that surface. Instead, we spec- 
ify the appropriate boundary conditions, thus linking the 
perturbation of the HI envelope to that of the filament. 



4.1 The Inner Boundary: 







The boundary condition at the radial centre of the filament 
(r = 0) is trivial for axisymmetric modes (m — 0). Mass 
conservation demands that poUr — > 0, and the perturbed 
gravitational field d^i/dr —> since the internal mass per 
unit length vanishes as r ^ 0. 



4.2 The Outer Boundary: r ^ 00 

The outer boundary conditions are that all components of 
the perturbation vanish at infinity; thus, all components of 
the momentum density povi, and the perturbed potential 
(fii, must vanish. 

It would be difficult to specify these boundary condi- 
tions on a uniform grid. Therefore, we employ a non-uniform 
radial grid spacing defined by the transformation 



Stan 5 



(14) 
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where 5 is a constant scale factor; thus radial infinity is 
mapped to ^ = 7r/2. In practice, we vary ^ over the range 
5i < ^ < 7r/2 — 82 for small 5i and 82, which define the in- 
ner and outer boundaries. Thus, our numerical grid extends 
fron nearly zero radius to essentially infinite radius. Our 
transformation has the benefit of excellent dynamic range; 
by choosing 5* appropriately, we can arrange to have most 
of the grid points (approximately 2/3) inside of the molecu- 
lar filament, where good resolution is most important, with 
progressively fewer points in the HI envelope as r ^ cxa. In 
practice, we find that the eigensystem is quite insensitive to 
the choice of ^2; there is no detectable change in the eigen- 
value — tj^ for 82 correspanding to outer radii greater than 
~ 100 ro, where ro is the core radius defined in Paper I for 
filamentary clouds. 

4.3 The Molecular Filament/HI Envelope 
Internal Boundary 

The most complex boundary in the problem occurs at the 
interface between the molecular and atomic gas. This surface 
is a contact discontinuity which moves freely, but through 
which no material may pass. No material may cross the 
boundary since this would involve a phase transition be- 
tween atomic and molecular gas, which requires much longer 
than the dynamical timescales relevant to our problem. The 
boundary conditions at the interface are written as follows: 



[$ 


] = 





(15) 


[V'I'-n 


] = 





(16) 


[B-n 


] = 





(17) 


[-(-+^)^+S(--)] 


= 


0, 


(18) 



where square brackets denote the jump in a quantity across 
the boundary. We note that these jump conditions apply in a 
Lagrangian frame that is co-moving with the deformed sur- 
face. Equations na and [19 demand that the gravitational po- 
tential and its first derivative (the gravitational field) must 
be continuous aross the interface. Equation nj\ is the usual 
condition on the normal magnetic field component from elec- 
tromagnetic theory; it is easily derived from the divergence- 
free condition of the magnetic field. The final condition, 
equation hq states that the normal component of the total 
stress is continuous across the boundary. 

Appendix H derives the explicit forms of these equa- 
tions in terms of the components of ^ (equation |l3| ), and 
Appendix O shows how the boundary conditions can be 
included in the matrix eigensystem given by equation h2l 
Equation bsl gives the final form of our eigensystem, which 
we solve, in Section H, for various equilibrium configurations. 
In general, equatioii|C5| has 4 x (A*' — 1) eigenvalues, where 
N is the size of each block matrix, but most of these rep- 
resent stable MHD waves. We generally solve only for the 
dominant mode, since fragmentation will be dominated by 
the fastest growing axisymmetric instability (with largest 
positive — oj^). For this reason, our problem is well suited 
for iterative methods. We solve equation fca using the well- 
known method of shifted inverse iteration (cf. Nakamura 
1991). 



5 DISPERSION RELATIONS AND 
EIGENMODES 

For any equilibrium state (which may be prepared using the 
formulation of Paper I) , we may determine a dispersion rela- 



tion for the dominant mode by solving equation C5 for 



as we vary the wave number k^. The wave number specifies 
the wavelength A over which the instability operates; 

27r 



A 



(19) 



We shall be primarily concerned with the following 
characteristics of our dispersion curves. 1) kz,max is the 
wavenumber of maximum instability, which may determine 
the scale for the separation of fragments by equation [19. 2) 
—Lu^a^ is the maximum squared growth rate for unstable 
modes, and determines the timescale on which the instabil- 
ity operated. 3) kz,crit is the maximum wave number for 
which a mode is unstable. Therefore, it sets the minimum 
length scale of the instability. We shall often write kz,max, 



and k. 



and 



-LU, 



T^ax in dimensionless form, defined by 



fOl^z^max 

2 



<^max (47rGpc 



(20) 



where tildes are reserved for dimensionless quantities for the 
remainder of this paper. 



6 TESTS OF THE NUMERICAL CODE 

We have tested our code by reproducing some of the dis- 
persion curves given by Nagasawa (1987) and Nakamura 
(1993). These test problems include the following: 1) the 
untruncated Ostriker solution (Nakamura), 2) untruncated 
filaments with purely poloidal, purely toroidal, and helical 
fields (Nakamura), and 3) pressure truncated filaments with 
constant poloidal fields (Nagasawa). 

Nagasawa assumes that the external medium is non- 
conducting and infinitely hot, with constant pressure and 
zero density. Our results are indistinguishable from his when 
we use his boundary conditions. When we include the efi'ects 
of perfect MHD (infinite conductivity) and finite density in 
the external medium, our results converge towards his in 
the limit of high velocity dispersion (Te ~ 10— 100. At lower, 
and more realistic, velocity dispersions {as « 5), we find 
that the instability is slightly decreased, although the gen- 
eral character of his curves is preserved. These effects are 
further discussed in Section m. 



7 STABILITY OF THE PRESSURE 

TRUNCATED OSTRIKER SOLUTION 

In this Section, we systematically examine the effects of the 
density and pressure of the external medium on the sta- 
bility of hydrostatic filaments. The equilibrium state is the 
isothermal Ostriker solution discussed in Paper I; 



P = 



[l + ir/roW 



(21) 



which may be truncated at any desired radius by external 
pressure (see Paper I). 
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Two independent effects are shown in Figure nl; curves 
1 to 3 sliow tlie effect of varying tlie density of the external 
medium, while the solid curves (curves 4 to 8) shows the 
effect of varying the external pressure, or equivalently the 
mass per unit length. 

We first consider the effect of varying the density of the 
HI envelope. Equation H gives the density, as a function of 
radius, for our polytropic envelope. We note that the density 
is controlled, mainly, by the velocity dispersion ae,s just out- 
side of the molecular filament. The dashed curves (1 to 3) in 
Figure hi have been computed using using the same equilib- 
rium solution for the molecular filament (see Table hi), but 
different values for ae,s- For all three curves, we have as- 
sumed that Ps/{P) = 0.801 and m/rriyir = 0.199. We find 
that envelopes with lower velocity dispersions and higher 
densities have a slightly stabilizing effect. Although we have 
taken the external medium to be non-self-gravitating, it re- 
sponds dynamically and self-consistently to the gravitational 
field and motions of the molecular filament. As the filament 
fragments, the deformation of the surface and the changes in 
the gravitational potential induce motions in the surround- 
ing gas. These motions slightly resist the fragmentation of 
the filament, because of the finite inertial mass of the en- 
velope. We have also computed one solution, shown as the 
curve numbered 1, where we have assumed that the exter- 
nal medium is infinitely hot, and of infinitely low density, 
which is consistent with Nagasawa's (1987) treatment of the 
boundary conditions (see Section H) . It is apparent that dis- 
persion relations computed with fully self-consistent bound- 
ary conditions converge to this curve in the limit of high 
velocity dispersion. In practice, we find essentially no differ- 
ence when (7e/(7c ~ 10. 

We now consider the effect of the external pressure, or 
equivalently, the effect of the mass per unit length. The virial 
equation from Paper I demonstrates that pressure trunca- 
tion is equivalent to a reduction of the mass per unit length 
for hydrostatic filaments; 

^ — M- '-' 

where m^ir = 2{a^)/G is the virial mass per unit length. 
Thus, Ps/{P) and m/ravir cannot be varied independently 
for unmagnetized filaments. The solid lines in Figure pj 
(curves 4 to 8) show dispersion relations for the truncated 
Ostriker solution, where we have assumed that the velocity 
dispersion of the HI envelope just outside of the filament is 
five times that of the molecular gas (see Table hi) . We find 
that more severely pressure truncated filaments, with higher 
Ps/{P) and lower m/rrivir, are less unstable than more ex- 
tended filaments. Thus, pressure truncation suppresses the 
gravitational instability of the filament. This result is easy to 
understand. Decreasing the mass per unit length decreases 
the gravitational force that drives the instability. Therefore, 
pressure dominated filaments, with external pressure compa- 
rable to the central pressure, are more stable than those that 
are dominated by self-gravity (eg. untruncated filaments). 

Figure H shows how the mass per unit length affects 
the wave number k^^rnax and growth rate —uj'^ax of the 
most unstable mode, as well as the critical wave number 
kz,crit- We find that the growth rate is significantly sup- 
pressed when m/rriTjir ~ 0.8. The wave numbers remain es- 
sentially constant until m/rrivir ~ 0.3. For smaller masses 
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Figure 1. We have computed dispersion curves for truncated, un- 
magnetized, isothermal filaments. The external medium has finite 
density and a velocity dispersion of Oe = 5(t for the solid curves. 
The curve shown as dot-dashed has an infinitely hot, zero density, 
external medium. The dashed curves have been computed using 
external media with velocity dispersions a^ = lOcr for the upper 
curve, and Ce = 2.5o" for the lower curve. Table hi summarizes the 
results shown here. 





m/m^ir 


Ps/{P) 


C 


f^z,max 


-'^Irax 


'^z,crit 


1 


0.199 


0.801 


0.472 


0.149 


0.063 


0.866 


] 


2 


0.199 


0.801 


0.47 


0.149 


0.0616 


0.861 J 


3 


0.199 


0.801 


0.432 


0.149 


0.0467 


0.791 




4 


1 





oo 


0.285 


0.117 


0.564 


5 


0.597 


0.403 


0.537 


0.313 


0.0882 


0.578 


6 


0.398 


0.602 


0.362 


0.35 


0.0687 


0.64 


7 


0.199 


0.801 


0.149 


0.462 


0.0581 


0.847 


8 


0.0995 


0.901 


-0.027 


0.639 


0.0549 


1.17 



" Infinitely hot external medium. 
^ "Hot external medium": <Je = lOcr 



CTe : 

"Cold external medium": Ce 



2.5a 



Table 1. We give the the wave numbers and growth rates for the 
fastest growing modes of truncated, unmagnetized filaments (See 
Figure ih. We also give the critical growth rate, beyond which the 
filament is stable against axisymmetric perturbations. The tildes 
indicate dimensionless variables, as defined by equations EOl The 
external medium has finite density and a velocity dispersion of 
cTg = 5(T, except where otherwise indicated in the footnotes. 



per unit length, the wave number increases, corresponding 
to fragmentation on smaller length scales. 

Two examples of eigenmodes are shown in figures H and 
W for unmagnetized filaments; the first shows the fragmenta- 
tion of the untruncated Ostriker solution, while the second 
is significantly truncated {m/mvir = 0.199). These figures 
show a single wavelength in the periodic fragmentation of a 
filament. We note that the size of the perturbation is greatly 
exaggerated in these figures, where we typically add a per- 
turbation of ~ 33% (in peak density) to the unperturbed 
solution. 
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m/m 



Figure 2. We plot the wave numbers k^^max and k^^crit-, as well 
as the squared growth rate —u)'^ of the most unstable mode, as a 
function of the vajm^i^ . 



The arrows in figures H and Ml show the velocity field, 
but the scaling is logarithmic due to the large range of ve- 
locities that must be shown. Specifically, arrow length oc 
logio{i + av/ac), where we choose a « 50. Both cases have 
purely poloidal velocity fields. The untruncated case shows 
an asymmetric infall onto the fragment, with almost purely 
radial infall near the fragments, and some radially outward 
motions between them. On the other hand, the truncated 
model shows infall along the axis, but the motions are actu- 
ally radially outward in the vicinity of the fragment. This is 
due simply to the outward bulging of the fragment, which 
can occur because there is no infalling material approaching 
from the radial direction. 

There are interesting motions in the gas outside of the 
truncated filament shown in Figure m We observe inward 
motions, that are driven by external pressure crushing the 
less dense parts of the filament between fragments. At the 
same time, the bulging of the filament near the fragments 
drives a radially outward flow in the surrounding gas. We 
emphasize that no mass is actually exchanged between the 
molecular filament and the surrounding atomic gas. Gas 
flows towards the fragments along the fllament axis, which 



N 




Figure 3. We show the eigenmode corresponding to the most 
unstable mode of the untruncated Ostriker solution (rofcz.maa; = 
0.285, -uj'^{4ttGpc)~'^ = 0.117). The contours show the density, 
while the arrows show the velocity field. We note that the segment 
shown is just one wavelength in the periodically fragmenting fila- 
ment. The arrows are scaled logarithmically with the velocity, as 
discussed in the text. 



then bulge outwards slightly. This induces circulation mo- 
tions in the surrounding envelope in which gas flows away 
from the fragments and towards the increasingly evacuated 
inter-fragment regions. 



8 STABILITY OF TRUNCATED FILAMENTS 
WITH PURELY POLOIDAL AND PURELY 
TOROIDAL FIELDS 

Following Paper I, we study the effects of poloidal and 
toroidal fields separately before considering the more gen- 
eral case of helical fields in Section bI We note that none of 
the equilibrium models used in this Section fall within the 
range of observationally allowed models from Paper I. These 
models are useful, however, in that they provide insight into 
the roles played by the field components in more general 
helical field models. 



8.1 Poloidal Field Results 

Figure M shows a sequence of dispersion curves, in which 
we have varied the poloidal fiux to mass ratio Fz while 
holding the mass per unit length constant at m = 5 mo 
(m/myir = 0.199). We observe that the poloidal field has 
a stabilizing effect on the filament, but that the stabiliza- 
tion saturates for V^ ~ 10. The critical and most unsta- 
ble wave numbers are decreased by the poloidal field, but 
this effect also saturates when F^ ~ 10. These results agree 
qualitatively with those obtained by Nagasawa (1987) for 
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Figure 4. We show the eigenmode corresponding to the most 
unstable mode of a truncated Ostriker model with mass per unit 
length m/mo = 5 {m/myir = 0.199). For this mode, rok^^max = 
0.462 and -uj^i'iwGpc)''^ = 0.449. 



pressure truncated isothermal filaments threaded by purely 
poloidal fields, although the degree of stabilization is greater 
in our model. Chandrasekhar and Fermi (1953) also find that 
poloidal fields stabilize incompressible filaments, but they do 
not find any saturation of the stabilizing effect. This is espe- 
cially apparent in Figure H, where we have plotted kz^max, 
kz,crit, and —lo'^^^ as a funtion of F^. 

Figures M and M show two eigenmodes corresponding, 
respectively, to a weak (F^ — 1) and relatively strong (F^ — 
5) poloidal field. Naturally, as the poloidal flux to mass ratio 
increases, and the field strengthens, the gas becomes more 
tightly constrained to move only along the field lines. When 
Fz ~ 10, the motions are almost entirely poloidal; further 
increasing r2 has no significant effect on the motions, and, 
hence, no significant effect on the dispersion curves. 




Figure 5. We have computed dispersion curves for truncated 
filaments with poloidal field only. In all cases, the mass per unit 
length is 5mo {m/myir = 0.199), which falls within the range 
of observationally allowed masses per unit length found in Paper 
I. We have also assumed an external medium which is perfectly 
conducting and of finite density, with cTe = 5cr. The labeling^ of 
the dispersion curves corresponds to the numbering in Table ti. 



C Ps/{P) 



^z.max 



~,2 



1 





0.149 


0.801 


0.462 


0.0581 


0.846 


2 


1 


0.146 


0.813 


0.421 


0.0513 


0.751 


3 


2.5 


0.133 


0.859 


0.325 


0.0348 


0.556 


4 


5 


0.116 


0.929 


0.254 


0.023 


0.428 


5 


10 


0.106 


0.977 


0.225 


0.0184 


0.377 


6 


25 


0.102 


0.996 


0.216 


0.017 


0.362 


7 


50 


0.101 


0.999 


0.214 


0.0168 


0.359 



Table 2. We give the the wave numbers and growth rates for the 
fastest growing modes of truncated filaments with poloidal field 
only (See Figure pi). We also give the critical growth rate, beyond 
which the filament is stable against axisymmetric perturbations. 
In all cases, the mass per unit length is 5mo {m/m^ir = 0.199), 
which falls within the range of observationally allowed masses per 
unit length found in Paper I. We have also assumed an external 
medium which is perfectly conducting and of finite density, with 
(Te = 5cr. 



8.2 Toroidal Field Results 

Finally we consider models in which the magnetic field is 
purely toroidal. We stress that these models are illustra- 
tive only, and certainly do not represent realistic models 
of filamentary clouds. Figure ^ shows a sequence of disper- 
sion curves, in which we have varied the toroidal fiux to 
mass ratio F^ while holding the mass per unit length con- 
stant at 7n = 5 mo {m/m-uir = 0.199). We observe that 
the toroidal field strongly stabilizes the cloud against frag- 
mentation. Contrary to the effects of the poloidal field, the 
stabilization does not saturate; in fact, the models are com- 
pletely stable against fragmentation when F^ ~ 10. Purely 
toroidal fields have a stabilizing effect because fragmenta- 
tion requires a flow of material along the axis of the flla- 
ment, which must result in substantial compression of the 
toroidal flux tubes. This results in a gradient in the magnetic 
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0.8 



T3 



= 0.6 



ro.4 





Figure 6. We plot the wave numbers k^^max and k^^crit^ as well 
as the squared growth rate —w^ of the most unstable mode, as a 
function of the F, . 



pressure I/Stt dB^/dz which resists this motion. When the 
toroidal flux to mass ratio F^ is sufHciently high, the toroidal 
flux tubes can resist compression and arrest fragmentation 
altogether. The effects of the varying F^ are most readily ap- 
parent in Figure |lC|, where we have plotted kz,max, kz,crit, 
and —LOmax a-s a function of of Fa. We have shown an ex- 
ample of an eigenmode in flgures hU where r0 — 5. In both 
cases, the toroidal field becomes strongest near the frag- 
ments. The gas motions are purely poloidal in both cases, 
with a "reverse flow" just outside of the filament. As in the 
purely hydrodynamic case, the gas is pushed radially out- 
ward near the fragment, and fiows toward the empty regions 
between them. However, the pinch of the toroidal field re- 
stricts the motions to a narrow band just outside the fila- 
ment. 





Figure 7. We show the eigenmode corresponding to the most 
unstable mode of a truncated magnetized model with a weak 
poloidal field (Fz = 1), and mass per unit length m/mo = 5 
(m/myir = 0.199). For this mode, rokz,max = 0.421 and 
—Ltj^{4:nGpc)~^ = 0.0512. The figures represent a) (left panel) 
density contours with superimposed poloidal velocity vectors. We 
note that there is no toroidal velocity when the field is purely 
poloidal. and b) (right panel) The magnetic field. We always rep- 
resent the magnetic field as a split-frame figure, with poloidal field 
vectors shown on the right, and contours to represent the toroidal 
field on the left. Since no toroidal field is generated in this case, 
the left side of the figure is left blank. 





r^ 
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Ps/{P) 


f^z,max 


-,2 


kz,crit 


1 





0.149 


0.801 


0.462 


0.0581 


0.846 


2 


1 


0.162 


0.764 


0.427 


0.0499 


0.787 


3 


2 


0.199 


0.683 


0.358 


0.0357 


0.655 


4 


3 


0.253 


0.598 


0.284 


0.023 


0.513 


5 


4 


0.318 


0.523 


0.218 


0.0139 


0.39 


6 


5 


0.392 


0.462 


0.165 


0.00811 


0.293 


7 


6 


0.471 


0.413 


0.123 


0.0046 


0.217 


8 


7 


0.555 


0.372 


0.0914 


0.00256 


0.16 



Table 3. We give the the wave numbers and growth rates for the 
fastest growing modes of truncated filaments with toroidal field 
only (See Figure W. We also give the critical growth rate, beyond 
which the filament is stable against axisymmetric perturbations. 
In all cases, the mass per unit length is 5mo {'m/'mvir = 0.199), 
which falls within the range of observationally allowed masses per 
unit length found in Paper I. We have also assumed an external 
medium which is perfectly conducting and of finite density, with 
(Je = 5cr. 
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Figure 8. We show the eigenmode corresponding to the most 
unstable mode of a truncated magnetized model with a rela- 
tively strong poloidal field (F^ = 5), and mass per unit length 
m/mo = 5 (m/ra^iT = 0.199). For this mode, rdk^^max = 0.225 
and -a)2(47rGpc)~^ = 0.0228. The figures represent a) (left 
panel) density contours with superimposed poloidal velocity vec- 
tors. We note that there is no toroidal velocity when the field 
is purely poloidal. and b) (right panel) The magnetic field. We 
always represent the magnetic field as a split-frame figure, with 
poloidal field vectors shown on the right, and contours to repre- 
sent the toroidal field on the left. Since no toroidal field is gener- 
ated in this case, the left side of the figure is left blank. 



9 STABILITY OF FILAMENTARY CLOUDS 
WITH HELICAL MAGNETIC FIELDS 

In this Section, we consider the stability of models that best 
agree with the observational constraints found in Paper I: 



0.11 < 



< 0.43 



0.012 < jpr < 0.75. 
0.2 < X < 5, 



where the reader is refered to Paper I for a description of 
the various quantities. We used a Monte Carlo exploration 
of our three dimensional (F^, F^, and C) parameter space 
to find a set of models that agree with these constraints. We 
use the same set of models here to determine the stability 
properties of filamentary clouds threaded by helical fields. It 
is not necessary to show full dispersion curves for all mod- 
els, since most of the useful information is contained in the 




Figure 9. We have computed dispersion curves for truncated 
filaments with toroidal field only. In all cases, the mass per unit 
length is 5mo {m/myir = 0.199), which falls within the range 
of observationally allowed masses per unit length found in Paper 
I. We have also assumed an external medium which is perfectly 
conducting and of finite density, with cTe = 5a. The labeline_of 
the dispersion curves corresponds to the numbering in Table a. 



dispersion parameters k^^max, kz^crit, and 

try to determine which combination of model parameters 

Fz, r^, and C controls the stability properties. 

We find that the dispersion parameters kz,max, k^^cru, 
and —uj^^^ are best correlated with the ratio F^/F^. Figures 
n3 to nil show scatter plots of these dispersion parameters as 
functions of V^/Tz- It is obvious from these figures that we 
have found two types of unstable mode in our calculation; 
the two ty pes of unstable mode are mostly well separated 
in Figures Il3 to 14, but join together when P^/Pz ~ 2. 
The first, which occurs when P^/Pz is small, is a gravity- 
driven instability, which is analogous to the modes found for 
purely hydrodynamic models (Section M), and MHD models 
in which the field is purely poloidal or toroidal (Section pf). 
We note that increasing P^/Pz substantially decreases the 
growth rate of the instability. This is in accord with the 
findings of Section N, where we found that increasing P^ and 
Fz both stabilize gravity driven modes, but F^ is much more 
effective. In fact, virtually all of the gravity-driven modes 
in Figure [13 have growth rates that are significantly lower 
than the growth rates of unmagnetized filaments, as well as 
filaments with purely poloidal fields (shown in figures hi and 
y respectively). Typically, we find 

~ 0.03 





0.0005 


^ 


^^max 


(23) 


0.1 


< 


l^z^max 


(24) 


0.2 


< 


i^z^crit 



-0.2, 

0.6, (25) 

where the tildes refers to the dimensionless growth rates and 
wave numbers, as defined in equat ion [20| . 

By taking the inverse of y/~u'^, we find a growth 
timescale of 



Tjr 



1.8 



0.01 



-1/2 



V 104 cm-3y " 



(26) 



where ric is the central number density of the filament. We 
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Figure 10. We plot the wave numbers k^^max 
as the squared growth rate 
function of the F^. 



and kz^crit^ as well 
of the most unstable mode, as a 



have chosen a fiducial central density of 10 cm~ because 
radially extended (0.5 pc diameter) filamentary structure is 
clearly visible in C^^O maps of Taurus (Onishi et al. 1998). 
C^^O molecules require a density of at least ~ 2 x 10^ cm'^^ 
for excitation, and our models predict that the central den- 
sity should be several times higher than the bulk of the fil- 
ament. Nevertheless, the central densities of filaments have 
not yet been accurately measured, so this fiducial central 
density should be treated with caution. The radial signal 
crossing time is approximately given by 



^ = 0.49 
a 



Rs 



0.25 pc 



0.5 km s 



Myr. 



(27) 



Since the signal crossing time is comparable to, and probably 
slightly longer than the fragmentation timescale, we expect 
filamentary clouds to achieve radial quasi-equilibrium before 
fragmentation destroys the filament in a few times Tfrag', 
thus, our analysis is consistent with the assumption of quasi- 
equilibrium used in Paper L 

The second type of unstable mode is driven by the 
magnetic field. These MHD-driven instabilities are triggered 
when r^/Fz ~ 2. Figures hd to hj show that this mode 
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Figure 11. We show the eigenmode corresponding to the most 
unstable mode of a truncated magnetized model with a toroidal 
field (r^ = 5), and mass per unit length m/mo = 5 (m/myir = 
0.199). For this mode, k^^„^ax = 0.165 and -uj'^{4.nGpc)^'^ = 
0.0081. The figures represent a) (left panel) density contours with 
superimposed poloidal velocity vectors. We note that there is no 
toroidal velocity when the field is purely toroidal, and b) (right 
panel) The magnetic field. We always represent the magnetic field 
as a split-frame figure, with poloidal field vectors shown on the 
right, and contours to represent the toroidal field on the left. 
Since no poloidal field is generated in this case, the right side of 
the figure is left blank. 



joins onto the gravity-driven modes, but extends to very 
high growth rates and very large wave numbers. Since large 
wave numbers correspond to small wavelengths, gravity can- 
not be important, as a driving force, for these modes. They 
are, however, most unstable when the toroidal fiux wrap- 
ping the filament is large, compared to the poloidal fiux. 
This suggests that the mode is an MHD instability. In fact 
the criterion that we have found is analogous to the famous 
stability criterion for MHD "sausage" modes: 



B^s > V2B, 



(28) 



(cf. Jackson 1975). The analogy is far from perfect because 
equation |2^ applies to the case of a uniform plasma cylin- 
der threaded by a constant poloidal magnetic field, and sur- 
rounded by a toroidal field driven by a thin axial current 
sheet along the surface of the plasma. However, the axial 
current distribution in our case is distributed throughout 
the filament. Nevertheless, the similarity between equation 
ps and our own instability criterion is very suggestive. 
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Figure 12. The most unstable wave numbers for a random sam- 
ple of models that agree with our constraints. 



Figure 14. The critical wave numbers for a random sample of 
models that agree with our constraints. 
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Figure 13. The growth rates for a random sample of models that 
agree with our constraints. 



We do not find MHD-driven instabilities for models 
with purely toroidal fields. We note that the stability cri- 
terion for a non-self-gravitating plasma column pinched by 
a purely toroidal field can be expressed analytically as 



dlnP 



d In r- 2 + 13 



(29) 



for an isothermal perturbation, where f} is the usual plasma 
P, defined by /3 = SnP/Bl (cf. Sitenko & Malnev 1995). We 
have checked, numerically, that our equilibria obey equa- 
tion p9l and, therefore, should be stable against "sausage" 
modes. It may appear that equation hs is violated for any 
model with a purely toroidal field, but we remind the reader 
that equation ES strictly applies to situations in which the 
toroidal field is generated by a current sheet at the plasma 
surface, and thus resides entirely outside of the plasma col- 
umn. Equation M is, however, appropriate for the more dis- 
tributed fields in our models. 

It seems paradoxical that models with purely toroidal 
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fields are stable against "sausage" modes, while those with 
a poloidal field component may be unstable. However, we 
note that helical fields are topologically different from purely 
toroidal fields. Flux tubes in purely toroidal models form 
closed loops, while the loops of toroidal fiux are all linked in 
the case of a helical field. When the field is helical, twisting 
motions in the filament may enhance the toroidal field, while 
no amount of twisting can amplify a purely toroidal field. 
Thus, we argue that it is at least possible for modes to exist 
that simultaneously generate twisting motions and toroidal 
field, which could destabilize the filament. 

We note that the gravity-driven modes join continu- 
ously with the MHD-driven modes in Figures hd to h4: 



intersection of the two branches, when V^/Tz 



the 
repre- 
sents a transition between gravity and the toroidally domi- 
nated magnetic field as the main driving force of the insta- 
bilty. Some of the modes that we have labeled as magnetic 
have growth rates that are similar to gravity modes, ex- 
cept that the wavelength is somewhat shorter. These may 
represent reasonable modes of fragmentation for filamen- 
tary clouds. However, the more unstable modes of the MHD 
branch probably cannot, since they grow too rapidly and 
would probably disrupt the filament on timescales that are 
much shorter than the timescale on which equilibrium can 
be established (See equations |2q and ETI). Therefore, a sub- 
set of the models allowed by the observational constraints 
of FPl are actually highly unstable. 

In Figure [19, we show an example of an unstable 
gravity-driven mode. As in the previously considered cases 
of hydrodynamic filaments (Section m and filaments with 
only one field component (Section 0), gravity drives a purely 
poloidal fiow of gas towards the fragments. However, a heli- 
cal field results in a toroidal component of the velocity field 
which alternates in sign from one side of the fragment to 
the other. The reason is that the helical field exerts oppo- 
sitely directed torques on the gas fiowing towards the frag- 
ment from opposite directions. The mode shown in Figure hq 
is driven by the MHD "sausage" instability. The dominant 
wavelength of the instability is much shorter than either the 
filament radius or the dominant wavelength of the gravity- 
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driven modes; therefore gravity is relatively unimportant 
compared to the magnetic stresses. The actual structure of 
the mode is quite similar to that of the gravity-driven modes 
except that the mode is confined to the most central parts 
of the filament, where the B^ gradient is steepest. This is 
consistent with our claim that this mode is a "sausage" in- 
stability, because "sausage" modes are driven by outwardly 
increasing toroidal fields with strong gradients. The filament 
is crushed where the toroidal field is strongest, which forces 
gas out along the axis and towards the fragments. As in 
the case of the gravity-driven modes, toroidal motions are 
generated as the helical field exerts torques on the gas. 

Figure O shows the effect of F^/Fz on the ratio of 
the fragmentation wavelength Xmax to the filament diam- 
eter D, which we define as twice the filament radius Rs; 
we note that gravity-driven and MHD-driven modes can- 
not be easily separated on this diagram. In all cases, we 
find that more toroidally dominated modes have smaller 
Xmax/D. The Schneider and Elmegreen (1979) Catalogue 
of Dark Globular Filaments shows that Xmax/D « 3.0 for 
most filaments in their sample. Interestingly, this value is 
consistent with F^/Fz ~ 2, which is near the transition from 
gravity to MHD-driven modes, and is near the lowest possi- 
ble growth rate. We postulate that most filamentary clouds 
may be observed to reside near this maximally stable con- 
figuration because such objects would survive the longest. 
Although suggestive, this conclusion will also depend on the 
non-axisymmetric stability of filamentary clouds, to be dis- 
cussed in a forthcoming paper. 

Finally, in Figure |l^, we show the masses of fragments 
that might form from filamentary clouds. We note that 
we have computed the fragment masses simply by multi- 
plying the equilibrium mass per unit length by the wave- 
length of the most unstable mode. This procedure implic- 
itly assumes that there is no sub- fragmentation as fragments 
evolve. Thus, these masses should be regarded with caution. 
We find that fragment masses, computed in this way, fall in 
a rather large range; in dimensionless units we find 



100 ~ Mfraai^nCypt'^a 



3/2 1/2 -3 < 



400 



(30) 



for gravity-driven modes, which have F^/Fz ^ 2. Thus, we 
obtain masses of 



Mfrag 



:41.5 



100 



0.5km 



ms-^J UO* cr 



-1/2 



10^ 



-10; 



10 



Mq. (31 



-10 



10 DISCUSSION AND SUMMARY 

We have determined how pressure truncation and helical 
fields affect the stability of our models against axisymmet- 
ric modes of fragmentation, which are ultimately responsible 
for the formation of cores. Our analysis difi'ers from previ- 
ous work on filaments in a number of respects. Firstly, and 
most importantly, we have made at least a first attempt 
to obervationally constrain our models; thus, our modes of 
fragmentation pertain to equilibria that are truncated by re- 
alistic external pressures and have magnetic fields that are 
consistent with the observations. Moreover, our models have 
density gradients of r~^* to r^'^ , in good agreement with 
observational data (Alves et al 1998, Lada et al. 1998). Sec- 
ondly, we have treated the external medium self-consistently 
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Figure 16. An example of an MHD-driven mode for a model 
with a helical magnetic field. For this mode, kz,max = 1-17 and 
—uj^{4:nGpc)~^ = 0.0390. The figures represent a) (top panel) 
density contours with superimposed poloidal velocity vectors, b) 
(middle panel) Contours of toroidal velocity. Light coloured con- 
tours represent positive v^, while dark coloured contours repre- 
sent negative 110. and c) (bottom panel) The magnetic field. We 
always represent the magnetic field as a split-frame figure, with 
poloidal field vectors shown on the right, and contours to repre- 
sent the toroidal field on the left. 
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Figure 15. An example of a gravity-driven mode for a model with a helical magnetic field. For this mode, k^^max = 0.199 and 
— a;^(47rGpc)~^ = 0.0081. The figures represent a) (left panel) density contours with superimposed poloidal velocity vectors, b) (middle 
panel) Contours of toroidal velocity. Light coloured contours represent positive v^, while dark coloured contours represent negative v^. 
and c) (right panel) The magnetic field. We always represent the magnetic field as a split-frame figure, with poloidal field vectors shown 
on the right, and contours to represent the toroidal field on the left. 



by considering it to be a perfectly conducting gas of finite 
density. 

We found, in Paper I, that filaments that are highly 
truncated by the external pressure have lower masses per 
unit length than untruncated filaments. We have shown, in 
Section M, that such pressure dominated equilibria are more 
stable than their less truncated gravitationally dominated 
counterparts. We also found that the finite density of the 
external medium, in our calculation, slightly stabilizes fila- 
ments due to inertial effects. This effect vanishes when the 
external medium is "hot", with a velocity dispersion ~ fO a. 

Both the poloidal and toroidal fields were shown to sta- 
bilize filaments against gravity-driven fragmentation in Sec- 



tion H. The poloidal field stabilizes the filament by resisting 
motions perpendicular to the poloidal field lines. This stabi- 
lization eventually saturates, when the field is strong enough 
{Vz ~ 10) to disallow radial motions altogether. Surpris- 
ingly, the toroidal field is much more effective at stabilizing 
filaments against gravitational fragmentation. As the fila- 
ments begins to fragment, the toroidal flux tubes are pushed 
together near the fragments; this results in a gradient in the 
magnetic pressure (I/Stt dBpf^^/dz) that opposes fragmen- 
tation. In the unreaUstic limit of purely toroidal fields, the 
filament is stable to fragmentation when F^ ~ 7. One might 
expect that these models would be highly unstable against 
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Figure 17. The ratio of fragment separation Xmax = 27r/fcmai 
to filament diameter D = 2Rs for a random sample of models 
that agree with our constraints. 



10' 



10 



10' 





■ ■;. 




' 




. ■■■■-•• -.V; 


■\«- 






■ '.'.ni?^' 


m- 




' 




% 

■■■4. 




' 






^ ■ ■ 




















10" 



10" 



r,/r. 



10' 



10' 



Figure 18. The expected mass of fragments for a random sample 

of models that agree with our constraints. Here we assume that 

It 

'-frag 



MHD "sausage" modes. However, we have verified that our 
equilibria satisfy the stability condition given in equation 

In Section HI, we demonstrated that our models with 
helical fields are subject to both gravity-driven and MHD- 
driven instabilities. The two types of modes blend together 
in Figures hd to M, with the transition from gravity-driven 
to MHD-driven largely determined by the ratio V^jV:^; when 
Tc^/Fz ~ 2, the instability is driven by gravity, with MHD- 
driven modes triggered only for more tightly wrapped helical 
fields. We provided evidence, in Section y, that the MHD- 
driven mode is probably is an MHD "sausage" mode. These 
instabilities can have extremely high growth rates compared 
to gravity driven modes; therefore, we may interpret our 
stability condition as an additional theoretical constraint on 
our models. We note that the most stable of all equilibria 
are near the transition from grivity-driven to MHD-driven 



modes. Thus, we find that there is a regime in which fil- 
amentary clouds fragment very slowly. For gravity-driven 
modes, we find that — 01^^^,, fez, max, and k^^crit fall within 
the ranges 



0.0005 
0.1 
0.2 



— UJ, 

rok. 
rokz 



2 

max 



{4ttGp. 
< 



0.03 



z.max 



0.3 



0.1 



(32) 



where ro is the core radius for filamentary clouds (see Paper 
I), and pc is the central density. We may also write the ex- 
pected wavelength for the separation of fragments and the 
growth timesacale as 



^max 



'''frag 



2.8 



1.8 



"02" 



) VO. 5 fern s-V UO^ cm-V 



-1/2 



0.01 



pc 



10-* 



-1/2 



Myr. (33) 



We note that the growth time for gravity-driven modes is 
much longer than the radial signal crossing time (equation 
P?!), on which radial equilibrium is established. Therefore, 
our stability analysis is consistent with the assumption of 
radial quasi-equilibrium used in Paper I. 

We find that the fragmentation of filamentary clouds 
with helical fields leads to a toroidal component of the ve- 
locity field, whose direction alternates with every half wave- 
length of the perturbation. These motions are comparable 
to, and in some cases exceed the poloidal velocities, at least 
in linear theory. These motions could, in principal, be de- 
tected, which would provide considerable evidence for helical 
fields. 

We find that the fiducial growth time for gravity-driven 
modes is on the order of 1.8 Myr, which is longer than the 
signal crossing time, on which radial equilibrium is estab- 
lished. Therefore, our stability analysis is consistent with 
the assumption of radial quasi-equilibrium used in Paper I. 

Throughout this analysis, we have focussed on axisym- 
metric modes that lead to fragmentation, ignoring any non- 
axisymmetric modes that might be present. Some of our 
models may be unstable to at least the m = 1 "kink" insta- 
bility. However, we do not expect gross instability for most 
models. While they do contain significant toroidal fields, the 
poloidal field is actually much stronger throughout most of 
the filament, particularly in the central regions, where B^ is 
maximal and B^ vanishes. This "backbone" of poloidal field 
should largely stabilize our models against kink modes. In 
any case, the non-axisymmetric modes of fragmentation will 
be addressed in a forthcoming paper. 
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APPENDIX A: THE EQUATIONS OF 
LINEARIZED MHD 

The equations of linearized, self-gravitating, perfect MHD 
that we solve are as follows: 
Momentum Equation: 

icjpovi + VPi +poV$i +piV$o 



— [(VxBo) X Bi + (VxBi) X Bo] = 
47r 



Mass conservation: 

iujpi + V ■ (povi) = 
Induction equation: 
itjBi = Vx(vixBo). 



Poisson's equation: 
Pi- 



V^$i 



Equation of state: 
Pi =7cropi, 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



where 7 is the adiabatic index of the gas (7 = 1 within the 
filament.) and o-q — Po/po is the one-dimensional velocity 
dispersion. 

By introducing a "modified" gravitational potential 
yji = ico^i, it is straightforward to combine equations W^ 
to AS to form the eigensysem given by equations M and fil 
Equations H and M represent the eigensystem, written en- 
tirely in terms of ipi and the components of pvi. We now 
define differential operators A, B, C, and I) in a way that 
is obvious from equations pi and H to write our eigensystem 
in closed symbolic form: 



LO poVl 





^(povi) + B(pi 
C(povi) + Dipi. 



(A6) 
(A7) 



We obtain an approximate matrix representation of op- 
erators A, B, C, and D by finite differencing over a one- 
dimensional grid of A'^ cells {N fa 500 usually). From this 
point forward, we shall assume that all operators have been 
finite differenced, and shall make no distinction between ma- 
trix and operator forms. We note that equation A7 applies 



to both the HI envelope and the molecular filament; thus, 
of the A'^ cells, some portion (usually ~ 2/3) are within the 
filament, while the remainder are in the HI envelope. Since 
the density is discontinuous at the interface, care must be 
taken not to difference any equations across the boundary. 
Section [J describes the boundary conditions that link the 
perturbation in the filament to the external medium. Matri- 
ces C and D are sparse tridiagonal N x N matrices, while 
A and -B each take the form of 3 x 3 blocks of N x N sparse 
tridiagonal sub-matrices (The 3x3 block structure occurs 
because these matrices operate on povi, which has 3 vec- 
tor components.). If we consider "eigenvector" ^ given by 
equation M, then equation A7 takes the form of the stan- 



dard eigenvalue problem given by equation h3, where 



L = 



A B 
C D 



(A8) 
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APPENDIX B: THE BOUNDARY CONDITIONS 
AT THE MOLECULAR FILAMENT/HI 
ENVELOPE INTERFACE 

Equations ha to hd give the boundary conditions at the in- 
terface between the molecular and atomic gas. We note that 
these jump conditions apply in a Lagrangian frame that is 
co-moving with the deformed surface. Equations |l5| and |l6| 
demand that the gravitational potential and its first deriva- 
tive (the gravitational field) must be continuous aross the 
interface. Equation |l7| is the usual condition on the normal 
magnetic field component from electromagnetic theory. It is 
easily derived from the divergence-free condition of the mag- 
netic field. The final condition, equation hq states that the 
normal component of the total stress is continuous across 
the boundary. 

We must now transform the boundary conditions (equa- 
tions hsl), into an Eulerian frame appropriate to our fixed 
grid. We assume that the deformed surface is defined by the 
equation 



r = Rs + ee 



i(ujt-\-m(i>-\-k^z) 



(Bl) 



where e << Rs for a small perturbation. We note that 
m = for the axisymmetric modes considered in this pa- 
per. However, we retain m in our equations, since we will 
examine non-axisymmetric modes in a forthcoming paper. 
To first order, the unit normal vector to this surface is given 
by 



Rs 



I + ikz 



i(^ijjt-\-m(ji-\-kzz) 



(B2) 



Since this is a contact discontinuity, and not a shock, the 
velocity field must be consistent with the motion of this 
surface; thus. 



dT ■ iii^t + mS+k^z) 

at 
which implies that 



Url 







(B3) 



(B4) 



in the Eulerian frame. Solving equation B3 for e and substi- 



tuting into equation B2, we obtain an expression for n that 



involves only the radial velocity: 



- Url 



m 7 fcz , 

-—-(j}+—Z 



(B5) 



We insert the expression for the normal vector n (equa- 
B5) into equations ha to hsl and evaluate all zeroth 



tion 

order quantities, at the position of the deformed surface, by 

a fir st order Taylor expansion. We also make use of equations 

ah 



to A 



in order to express all quantities in terms of the 
momentum density ip and the potential ip. After some alge- 
bra, we obtain the explicit form of the boundary conditions 
that apply at the surface of the filament: 



[f^\ = 



dr 



+ pOVr 



(P0«r 

PO 



= 0. 



-70-0 V- (povi) 



Bo _ / Bo 

-— ■ Vx — xpovi 
4vr V po 



(B6) 
(B7) 

(B8) 
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+ 



1 dPtot 
Po dr 



(povi] 



= 0. 



(B9) 
(BIO) 



For completeness, we also study the stability of un- 
truncated filaments, for which there is no external pressure. 
When untruncated filaments extend radially to infinity, we 
simply do not include an internal boundary. We also consider 
the limit of an infinitely hot, zero density external medium, 
fn these cases, we assume that the external medium is a 
non-conducting vacuum and use the boundary conditions 
prescribed by Nagasawa (f987). 



APPENDIX C: MATRIX REPRESENTATION 
OF THE BOUNDARY CONDITIONS 



The four boundary conditions given by equation B6 to BE 
must now be inserted into the eigensystem given by equa- 
tion hi Assuming that there are A'' elements in our compu- 
tational grid, with the interface between the molecular and 
atomic gas between elements I* and /* -I- 1, we may write 
a matrix representation of our boundary conditions in the 
form 



(CI) 



a** -1- 6* = 0, 


where 


" *i 


** = 


*7* 




*iV 



(C2) 



contains all of the boundary cells, and '^ now contains the 
remaining 4(A^ — 1) elements. We likewise separate our eigen- 
system (equation hd) into regular and boundary grid cells: 

(C3) 



-LJ ^ Li* -I-L2*. 

We may solve equation 
mentary matrix algebra: 



CI 



for the boundary values by ele- 



(C4) 



Substituting into equation C3, we eliminate 4 equations to 



obtain our eigensystem (of size 4{N — 1) x 4(A'' — f)), which 
takes into account all boundary conditions: 



(L2 



LlQ.^^b)*. 



(C5) 



This is the final form of our eigensystem, which we solve in 
Section 0. 
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